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', Abstract 

As Cosmic Microwave Background (CMB) measurements are becoming more am- 
bitious, the issue of foreground contamination is becoming more pressing. This is 
, especially true at the level of sensitivity, angular resolution and for the sky coverage 

of the planned space experiments MAP and Planck. We present in this paper 
an indicator of the accuracy of the separation of the CMB anisotropies from those 
induced by foregrounds. 

Of course, the outcome will depend on the spectral and spatial characteristics of 
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the sources of anisotropies. We thus start by summarising present knowledge on the 
spectral and spatial properties of Galactic foregrounds, point sources, and clusters 
of galaxies. This information comes in support of a modelling of the microwave sky 
including the relevant components. The accuracy indicator we introduce is based on 
a generalisation of the Wiener filtering method to multi-frequency, multi-resolution 
data. While the development and use of this indicator was prompted by the prepa- 
ration of the scientific case for the Planck satellite, it has broader application 
since it allows assessing the effective capabilities of an instrumental set-up once 
foregrounds are fully accounted for, with a view to enabling comparisons between 
different experimental arrangements. 

The real sky might well be different from the one assumed here, and the analysis 
method might not be in the end Wiener filtering, but this work still allow meaningful 
comparative studies. As a matter of examples, we compare the CMB reconstruction 
errors for the MAP and Planck space missions, as well as the robustness of the 
Planck outcome to possible failures of specific spectral channels or global variations 
of the detectors noise level across spectral channels. 
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1 Introduction 



Future Cosmic Microwave Background (hereafter CMB) anisotropy experi- 
ments should tightly constrain the parameters of any theory for the origin 
of those anisotropies. In the case of the Planck space mission of ESA (for- 
merly known as COBRAS/SAMBA), for inflation generated fluctuations, 
the value of the cosmological parameters should be constrained at the percent 
level [ Bond et al., 1997| , Zaldarriaga et ai., 1997, Etstathiou and Bond, 1998 



These estimates generally rely on assuming that the measurements at some 
frequencies are dominated by the primary cosmological signal, the other being 
dominated by foregrounds emissions and used to "clean" the former mea- 
surements from any substantial remaining contamination. Which frequency 
measurements to consider appears as a matter of art. 



In this paper, we propose a scheme for assessing quantitatively the effect of 
foregrounds from an assumed model of the microwave sky, a description of the 
experiment, and by supposing Wiener filtering is used to analyse the data. 
While linearly optimal, this separation method has also difficulties which we 
shall discuss later. Other non-linear methods are likely to be better suited 
for recovering the non-Gaussian spatial distributions of the foregrounds (and 
possibly of the CMB itself). But the linear nature of Wiener filtering allows 
straightforward calculations leading in particular to an effective spatial win- 
dow (or beam profile) and an effective noise level for the experiment, once 
foregrounds have been accounted for. This allows putting previous cosmolog- 
ical parameters accuracy forecasts on a firmer ground, and it is a useful tool 
for optimising experimental set-ups. 

Undoubtly new aspects of the foreground emission will be discovered in future 
experiments, and in particular by MAP and Planck themselves, but we hope 
that the present observations and understanding of foregrounds is already 
sufficient to allow building a plausible model of the microwave sky. The more 
so that we shall only need a rough statistical description to test the ability of 
planned experiments to disentangle foregrounds from true CMB anisotropies. 

This paper has three main parts. The first one is concerned with deriving a 
plausible microwave sky model adapted to our present purposes. The next 
one describes the Wiener filtering method in the context of CMB experiments 
and leads to the definition of a "quality factor" which assesses the effect of 
foreground contamination on the outcome of an experiment. The last part 
offers practical examples in the case of MAP and Planck by making use of 
our sky model and analysis method. More precisely, section 2 discusses the 
Galactic foregrounds, while section 3 discusses the extragalactic foregrounds 
due to infrared galaxies, radio-sources, and the Sunyaev-Zeldovich effect on 
clusters of galaxies. These various contributions are put into perspective with 
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respect to CMB anisotropies in § 4. Section 5 discusses component separation 
methods, Wiener filtering, and derives the properties of the proposed quality 
factor. Practical applications are considered in § 6. 

Before we go on, let us first establish notations. The amplitude of any scalar 
field (in particular the temperature anisotropy pattern) on the sphere, A(9, 4>), 
can be decomposed into spherical harmonics as 

A(6,<f ) ) = J2ae m Ye m (0A)- (1) 

The multipole moments, ae m , are independent for a statistically isotropic field, 
i.e. (aj m ) = CgSur S mm >; Ci is called the angular power spectrum and it is 
related to the variance of the distribution by 

c ^ {A i) = Yi *A±±C, (2) 
The angular power spectrum characterises completely a Gaussian field. 



2 Galactic Foregrounds 



The Galactic emissions are associated with dust, free-free emission from ionised 
gas, and synchrotron emission from relativistic electrons. We are still far from 
an unified model of the galactic emissions. On the other hand, the primary 
targets for CMB determinations are regions of weakest galactic emission, while 
stronger emitting regions are rather used for studies of the interstellar medium. 
The approach taken here is thus to find a model appropriate for the best half 
of the sky from the point of view of a CMB experiment. This is based on work 
done during the preparation of the scientific case for Planck (see Pouchet 



et al. 1997|, Bersanelli et al. 1996 1, AAQ 1998 



2.1 Dust emission 



This restriction to the best half of the sky results in considerable simplification 
for modelling the dust emission, since there is now converging evidence that 
the dust emission spectrum from high latitude regions with low HI column 
densities can be well approximated with a single dust temperature and v 2 
emissivity, with no clear evidence for a very cold dust component in these 
regions. 

The spectrum of the dust emission at millimetre and sub-millimetre wave- 
lengths has been measured by the Far-Infrared Absolute Spectrophotometer 
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(FIRAS) aboard the Cosmic Background Explorer (COBE) with a 7 
degree beam. In addition, several balloon experiments have detected the dust 
emission at high galactic latitude in a few photometric bands with angular 
resolution between 30 arcminutes and 1 degree (see for example [Fischer et al 



1995|| and |De Bernardis et al. |1992|| ). The spectrum of the Galaxy as a whole 



Wright et al., 199l| , |Reach et al., 1995H cannot be fitted with a single dust tem- 



perature and emissivity law. These authors have proposed a two temperatures 
fit including a very cold component with T ~ 7K. The interpretation of this 
result for the Galactic plane is not straightforward because dust along the line 
of sight is expected to spread over a rather wide range of temperatures just 
from the fact that the stellar radiation field varies widely from massive star 
forming regions to shielded regions in opaque molecular clouds. On the other 
hand, a consistent picture for the high galactic latitude regions has recently 
emerged. 



Boulanger et al. [1996J analysed the far-IR & sub-mm dust spectrum by se- 
lecting the fraction of the FIRAS map, north of —30° declination, where the 
H I emission is weaker than 250 K.km/ s (36% of the sky). The declination limit 
was set by the extent of the HI survey made with the Dwingeloo telescope 
Purton and Hartmann, 1994j| . For an optically thin emission this threshold 



corresponds to N(HI) = 5 x 10 20 cm~ 2 . Below this threshold, the correlation 
between the FIRAS and Dwingeloo values are tight. For larger HI column 
density the slope and scatter increase, which is probably related to the contri- 
bution of dust associated with molecular hydrogen since this column density 
threshold coincides with that inferred from UV absorption data for the pres- 
ence of H2 along the line of sight fSavage et al., 1977H . Poulanger et al. |1996|| 



found that the spectrum of the H I-correlated dust emission seen by FIRAS 
is well fitted by one single dust component with a temperature T = 17.5 K 
and an optical depth r/N H = 1.0 x 10- 25 (A/250 / um)" 



~ 2 cm 2 . 



The residuals to this one temperature fit already allow to put an upper limit on 
the sub-mm emission from a very cold component almost one order magnitude 
lower than the value claimed by [Reach et al. |1995|] . Assuming a temperature 
of 7 K for this purported component one can set an upper limit on the optical 
depth ratio r(7K)/r(17.5K) ~ 1. Pwek et al. |1997|| have also derived the 
emission spectrum for dust at Galactic latitudes larger than 40°, by using 
the spatial correlation of the FIRAS data with the 100/zm all-sky map from 
the Diffuse Infrared Background Experiment (DIRBE). They compared this 
spectrum with a dust model not including any cold component. The residuals 
of this comparison shows a small excess which could correspond to emission 
from a very cold component at a level just below the upper limit set by the 
Boulanger et al. |1996|| analysis. Both studies thus raise the question of the 
nature of the very cold component measured by [Reach et al. [1995 . 



Puget et al. [1996J found that the sub-mm residual to the one temperature fit 
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of Poulanger et al. |1996|| is isotropic over the sky. To be Galactic, the residual 
emission would have to originate in a halo large enough (> 50 kpc) not to con- 
tradict the observed lack of Galactic latitude and longitude dependence. Since 
such halos are not observed in external galaxies, Puget et al. |1996f suggested 
that the excess is the Extragalactic Background Light due to the integrated 
light of distant IR sources. Since then this cosmological interpretation has re- 
ceived additional support. First, the main uncertainty in these analysis was 
the difficulty to establish the contribution from dust in the ionised gas. |La- 
gachc [1998| j has extracted a galactic component likely to be associated with 
the ionised gas. Second, the analysis was redone QGuiderdoni et al. |1997|| , |La- 
gache 1 1998|1 , ILagache et al. |1998a| ) in a smaller fraction of the sky where the 
HI column density is even smaller (Nh i < 1 x 10 20 Hcm -2 ). In these regions, 
the H I-correlated dust emission is essentially negligible and the "residual" is 
actually the dominant component. This leads to a much cleaner determina- 
tion of the background spectrum which is slightly stronger than the original 
determination. As we shall see below in § 3.1, the cosmological interpretation 
fits well with the results of recent IR searches for the objects that cause this 
background. Finally, the recent analyses by [Schlegel et al. [1998| | and |Hauser 
et al. | 1998|| also detected a residual significantly higher than the one obtained 



by |Guiderdoni et al. |1997|| , |Lagache |1998|| , |Fixsen et al. 1 19981 and [Lagache 



et al. |1998a|| . The determination of the extragalactic component in (Schlegel 
et al. 1 1998|| is not very accurate from their own estimate. The difference with 
Hauser et al. |199"8|1 is mostly due to the fact that these authors do not find 



any evidence for the emission associated with ionised gas. 



In summary, there is now converging evidence that the dust emission spectrum 
from most of the high latitude regions with low H I column densities can be 
well approximated with a single dust temperature and v 2 emissivity with no 
evidence for a very cold dust component. The sub- mm excess seen in the 
FIRAS data is rather of extragalactic origin. Of course, this only holds strictly 
for the emissions averaged over the rather large lobe of the FIRAS instrument, 
and for regions of low column density. In denser areas like the Taurus molecular 
complex, there is a dust component seen by IRAS at 100/im and not at 60yum. 
This cold IRAS component is well correlated with dense molecular gas as 
traced by 13 CO emission [[Laureijs et al., 199 1| , [Abergel et al., 1994] , [Lagache 



et al., 1998b| ] . In the following, we assume that the CMB data analysis will be 
performed only in the simpler but large connected region of low HI column 
density. 



Concerning the scale dependence of the amplitude of the fluctuations, |Gautier 
et al. "1992 1 found that the power spectrum of the 100 /im IRAS data is 
decreasing approximately like the 3 rd power of the spatial frequency, £, down 



to the IRAS resolution of 4 arc minutes. More recently, |Wright |1998|| analysed 
the DIRBE data by two methods, and concluded that both give consistent 
results, with a high latitude dust emission spectrum, C(£), also cx £~ 3 in the 
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range 2 < t < 300, i.e. down to ^ ~ 60/ 1 ~ 0.2 deg (with C(IO) 1 / 2 ~ 
.22MJy/sr at |6| > 30). 



2.2 Free- free emission 



Observations of Hq, emission at high Galactic latitude and dispersion mea- 
surements in the direction of pulsars indicate that low density ionised gas (the 
Warm Ionised Medium, hereafter the WIM) accounts for about 25% of the gas 
in the Solar Neighbourhood fReynolds, 1989|| . The column density from the 



mid-plane is estimated to be in the range 0.8 to 1.4 x 10 cm" . Until recently, 
little was known about the spatial distribution of this gas but numerous Hq, 
observing programs are currently in progress. An important project is the 
northern sky survey started by Reynolds (WHAM Survey). This survey con- 
sist of H Q spectra obtained with a Fabry-Perot through a 1° aperture. The 
spectra will cover a radial velocity interval of 200 km/s centred near the LSR 
with a spectral resolution of 12 km/s and a sensitivity of 1 cm -6 pc (5a). 
Several other groups are conducting H a observations with wide field camera 
(10°) equipped with a CCD and a filter ||Gaustad et al., 1998|| . These surveys 
with an angular resolution of a few arc minutes should be quite complemen- 
tary to the WHAM data and should allow covering about 90 % of the sky. 
From these H Q observations one can directly estimate the free-free emission 
from the WIM (since both scale in the same way with the emission measure 
oc J nidi, see |Valls-Gabaud 1 19981 for a careful discussion). 



Kogut et al. [1996J have found a correlation at high latitudes between the 
DMR emission after subtraction of the CMB dipole and quadrupole and the 
DIRBE 240/im map. The amplitude of the correlated signal is plotted in 
Figure 3. The observed change of slope of the correlation coefficients at 90 
GHz is just what is expected from the contribution of the free-free emission 
as predicted by Pennett et al. |1994|| . The spectrum of this emission may be 



described by I v oc v 6 (the index is the best fit value of [Kogut et al. [19T?5|]). 



More recently Pe Oliveira-Costa et al. [1997] cross-correlated the Saskatoon 



data jTegmark et al., 1997|1 with the DIRBE data and also found a correlated 
component, with a normalisation in agreement with that of [Kogut et al. [1996 



This result thus indicates that the correlation found by [Kogut et al. 
persists at smaller angular scales. 



T997j 



Veeraraghavan and Davies [1997J used Ha maps of the North Celestial Pole 
(hereafter NCP) to determine more directly the spatial distribution of free-free 
emission on sub-degree scale. Their best fit estimate is = 1.3^oi7 £~ 2 - 27±om 
/iK 2 at 53 GHz, if they assume a gas temperature ~ 10 4 K. While this spec- 
trum is significantly flatter than the £~ 3 dust spectrum, the normalisation is 
also considerably lower than that deduced from [Kogut et al. [199B|. Indeed, 
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their predicted power at t = 300 is then a factor of 60 below that of the 
CO BE extrapolation if one assumes an £~ 3 spectrum for the free-free emis- 
sion. Additionally [Leitch et al. |1997f | found a strong correlation between their 



observations at 14.5 and 32 GHz towards the NCP and IRAS 100/xm emis- 
sion in the same field. However, starting from the corresponding H Q map, they 
discovered that this correlated emission was much too strong to be accounted 
for by the free-free emission of ~ 10 4 K gas. These and other results suggest 
that only part of the microwave emission which is correlated with the dust 
emission traced by DIRBE at 240/im may be attributed to free-free emission 
as traced by Ha. The missing part could be attributed to free- free emission 
uncorrelated with the dominating component at 240/im. 



Draine and Lazarian [1998J recently proposed a different interpretation for 
the observed correlation. They propose that the correlated component is pro- 
duced by electric dipole rotational emission from very small dust grains under 
normal interstellar conditions (the thermal emission at higher frequency com- 
ing from large grains). Given the probable discrepancy noted above between 
the level of the correlated component detected by [Kogut et al. 1 19961 an d the 



lower level of the free-free emission traced by Ha, it appears quite reasonable 
that a substantial fraction of this correlated component indeed comes from 
spinning dust grains. The spectral difference between free-free and spinning 
dust grains shows up at frequencies £ 30 GHz, i.e. outside of the range 
probed by MAP and Planck. They are thus spectrally equivalent for the 
2 experiments. However the spatial properties may be different. Indeed the 
spinning dust grains should be well mixed with the other grains resulting in 
tight correlations persisting at all spatial scales. The free-free emission might 
instead decorrelate from the dust emission at very small scales (since H II re- 
gions might be mostly distributed in "skins" surrounding dense H I clouds as 
proposed by |MCKee and Ostriker [1977] , a model supported by comparisons 



of IRAS and H a maps). 

In the following we shall assume for simplicity that in the best half of the 
sky (therefore avoiding dense HI clouds) the correlated microwave emission 
with I u oc z/~ - 16 in the MAP and Planck range is well traced at all relevant 
scales by the dust emission which dominates at high frequency and is well 
traced by HI. Still there may well be an additional free-free emission which 
is not well traced by the H I-correlated dust emission. Indeed the correlated 
emission detected by [Kogut et al. 1 19961 might be only part of the total signal 



with such a spectral signature that DMR detected. While all of this signal 
may be accounted for by the correlated component, the error bars are large 
enough that about half of the total signal might come from an uncorrelated 
component^. This is what we shall assume below. 



1 At 53 GHz, on the 10° scale, the total free-free like signal is AT// = 5.2±4.2^K, 
while the correlated signal is AT cor = 6.8±1.6//K. If we assume AT cor = 5.2/uK and 
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2.3 Synchrotron emission 



Away from the Galactic plane region, synchrotron emission is the dominant 
signal at frequencies below ~5 GHz, and it has become standard practice to 
use the low frequency surveys of |Haslam et al. |1982|| at 408 MHz and [Re 



|ich and Reich |1988| at 1420 MHz to estimate by extrapolation the level of 
synchrotron emission at the higher CMB frequencies. This technique is compli- 
cated by a number of factors. The synchrotron spectral index varies spatially 
due to the varying magnetic field strength in the Galaxy | Lawson et al., 19"%T| . 
It also steepens with frequency due to the increasing energy losses of the 
electrons. Although the former can be accounted for by deducing a spatially 
variable index from a comparison of the temperature at each point in the two 
low frequency surveys, there is no satisfactory information on the steepening 
of the spectrum at higher frequencies. As detailed by |Davies et al. 1 1996b 



techniques that involve using the 408 and 1420 MHz maps are subject to many 
uncertainties, including errors in the zero levels of the surveys, scanning errors 
in the maps, residual point sources and the basic difficulty of requiring a very 
large spectral extrapolation (over a decade in frequency) to characterise use- 
ful CMB observing frequencies. Moreover, the spatial information is limited 
by the finite resolution of the surveys: 0.85° FWHM in the case of the 408 
MHz map and 0.6° FWHM in the case of the 1420 MHz one. But additional 
information is available from existing CMB observing programs. 



For instance the Tenerife CMB project [[Hancock et al., 1994| , Pavies et al. 



1996a|| has observed a large fraction of the northern sky with high sensitivity at 
10.4 and 14.9 GHz (and also at 33 GHz in a smaller region). A joint likelihood 
analysis of all of the sky area implies a residual rms signal of 24/iK at 10 GHz 
and 20/xK at 15 GHz. Assuming that these best fit values are correct, one can 
derive spectral indices of a = —1.4 between 1.4 and 10.4 GHz and a = —1.0 
between 1.4 and 14.9 GHz. These values, which apply on scales of order 5°, 
are in agreement with those obtained from other observations at 5 GHz on 
~ 2° scales ||Jones et al., 1998" . 



At higher frequency, the lack of detectable cross-correlation between the Haslam 
data and the DMR data leads [Kogut et al. |1995fl to impose an upper limit 
of a = —0.9 for any extrapolation of the Haslam data in the millimetre wave- 
length range at scales larger than ~ 7°. In view of the other constraints at 
lower frequencies, it seems reasonable to assume that this spectral behaviour 
also holds at smaller scales. 



The spatial power spectrum of the synchrotron emission is not well known and 
despite the problems associated with the 408 and 1420 MHz maps it is best 
estimated from these. We have computed the power spectrum of the 1420 MHz 

AT uncor = AT cor , we have ATfj = 7.3/iK which is well within the allowed range. 
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map for the sky region discussed above. The results show that at £ ^ 100 the 
power spectrum falls off approximately as £~ 3 (i.e. with the same behaviour 
as the dust emission). 



2.4 A simple galactic model 



The galactic emission is strongly concentrated toward the galactic plane. But 
what is the geography of the galactic fluctuations? Do we expect large con- 
nected patches with low levels of fluctuations? What is the amplitude of fluc- 
tuations typical of the best half of the sky? 

To answer these questions, we created a galactic model valid at scale £ 1° from 
spectral extrapolations of spatial templates taken from existing observations. 
The 408 MHz full-sky map of Haslam is our template for the synchrotron 
emission, extrapolated to other frequencies with a spectral index of a = —0.9. 

We use the DIRBE 240/iin map as a template for the H I-correlated dust and 
the associated free-free emission detected by [Kogut et al. [1996] . As mentioned 



earlier, the global free- free like emission show that the HI correlated emission 
accounts for most (at least 50%), and maybe all, of the global signal. We con- 
servatively assume that there may be a second, H I-uncorrelated, component 
accounting for 5% of the total dust and 50% of the free-free emission. There 
is of course no known template for such a component, but we assume that it 
should have the same "texture" than the correlated component: we simulate it 
by using again the DIRBE 240 jum map, but North/South inverted (in galac- 
tic coordinates). While arbitrary, this choice preserves the expected latitude 
dependence of the emission, and the expected angular scale dependence of the 
fluctuations. The dust spectral behaviour is modelled as a single temperature 
component with = 18K and v 2 emissivity. We assume that the global free- 
free like emission will behave according to I v oc z/~ ' 16 , and normalise it to give 
a total rms fluctuation level of 7.3 /xK at 53 GHz. Since the weakest emissions 
in DIRBE are at the 1 MJy/sr level (and at the 10K level for Haslam), we 
immediately see that the corresponding fluctuations should be at the level of 
a few fiK around 100 GHz. 

Using this model, we can now create full sky maps of the total galactic emis- 
sion converted to an equivalent thermodynamical temperature at any chosen 
frequency. Note that at this resolution, extragalactic point sources and the SZ 
effect should contribute negligible signal. In order to obtain local estimates of 
the level of fluctuations, we have computed their rms amplitude (or), over 
square patches of 3 degrees on each side (i.e. containing about 11 beams of 
1 degree FWHM). Thus, this estimator only retains perturbations at angular 
scales roughly between 1 and 3 degrees, equivalent to restricting the contribu- 
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tion to the variance to a range of angular modes, with £ between ~ 20 and 60. 
Figure 1 shows an example of such a map at A = 2.4 mm (125 GHz). The two 
darkest shades of grey (a? ~ 10 /iK) delimit the area where galactic fluctu- 
ations would be comparable to those of a COBE-normalised CDM modeQ; 
they only represent 20% of the total sky area at this scale. This map pro- 
vides direct graphical evidence that a large fraction of the sky is quite "clean" 
around the degree scale ||Bouchet et al., 1996 b|,p|] . 



maps available in the full version at \f'tp : // ftp.iap.fr/pub/ fromjusers/ bouchet/ wiener7.ps. g^ 



Fig. 1. Top: map of the temperature fluctuations due to the galactic foregrounds 
at A = 2.4 mm. The fluctuation levels, denoted by 5 shades of grey, are obtained 
by computing the local variance over square patches with 3 degrees on a side (see 
text). From the lightest (small ctt) to the darkest regions close to the galactic plane, 
the levels correspond to ut < 1.0, 3.16, 10.0, 31.1//K, and <tt > 31.1/uK. As one can 
see on a map produced with the same method on a CMB fluctuations maps (in a 
standard CDM model), one would have ot ~ 10 uK, i.e. the galactic fluctuations 
should be at least 10 times smaller than the CMB ones (in that theory) in ~ 58% 
of the sky (at 2.4 mm). 



2 The parameters of this model are fi& = 0.05 for the baryonic abundance, £l c = 0.95 
for that of CDM, A = 0, Hq = 50 Km/s/Mpc, no reionisation, and only scalar modes 
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Cumulative Distribution of the Variance at 8.41 mm g D „„ = -2.00, T Eust = 18.00K. a F „ c = -0.15, a s , nc = -0.90 




Log a in K, tX Du3t = -2.00, V sl = 18.00K, o> re? = -0.10, a 3 = -0.90 Central v [GHz] of the bonds 

u t computed with a 60' FWHM beam on patches 174,906' on a side o.,. computed with a 60' FWHM beam on patches 174.906' on a side 



Fig. 2. a) Cumulative distribution functions of the galactic foreground fluctuations 
at A = 2.4 mm (125 GHz) in 3 degree patches for a 1 degree FWHM Gaussian 
beam. It gives directly the the fraction of the sky in which the rms fluctuation level 
is lower than a certain value. The CDM curves gives for comparison the expected 
range of fluctuations for a COBE normalised CDM model, b) Median fluctuation 
values (i.e. corresponding to 50% of the sky) for the galactic components around 
the degree scale as a function of frequency. The galactic fluctuations have a clear 
minimum around 100 GHz. 

Analysing figure 1 more quantitatively, figure 2. a shows the fraction of the sky 
for which the rms fluctuations of each of the foregrounds and of their total 
contribution are less than a given value at A = 2.4 mm. These curves are 
steep for low fluctuation levels; the best half of the sky (median value) is only 
a factor of ~ 4 noisier than the best regions, but a factor ~ 10 below the CDM 
level. The high level "plateaus" correspond to low galactic latitude regions. 
Indeed the maps tell us that that the "clean" sky essentially corresponds to 
the areas at high galactic latitudes, excepting only a few hot spots like the 
Magellanic clouds. Figure 2.b, which plots the median fluctuation values as a 
function of frequency, shows a clear minimum around 100 GHz, at a level lower 
than 10 % of the expected level of the CMB anisotropies. Thus measurements 
restricted to that frequency would already be sufficient to achieve better than 
10% accuracy around the degree scale (provided a few high latitude "hot or 
cold spots" are flagged out with higher and lower frequency channels). 

As was shown above, the angular power spectra C(l) of the galactic com- 
ponents all decrease strongly with £, approximately as C(£) oc £~ z . Smaller 
angular scales thus bring increasingly small contributions per logarithmic in- 
terval of £ to the variance, £{£ + 1)C(£): the galactic sky get smoother on 
smaller angular scales. We set the normalisation constants of the spectra to 
obtain a level of fluctuations representative of the best half of the sky. In prac- 
tice, this was achieved by generating an artificial map with C(£) = £~ 3 . We 
then measured its a? using our previous procedure of measuring the variance 
in 3 degree patches after convolving with a one degree beam. This gives the 
required normalisation factor of the power spectra by comparisons with the 
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Expected contributions in the model 




Contributions excluding the dipole and quadrupole 

Fig. 3. Standard deviation of the fluctuations in the synchrotron, H Icorrelated and 
uncorrelated emissions in our galactic model, as compared to the CMB fluctua- 
tions whose level was set by the DMR measurement, and the data points for the 
DMR-DIRBE correlated component found by Kogut et al. 95 (see text). A Gaussian 
beam of 7 degree FWHM was used. 

median values plotted in fig.2.b. We find that 

EC 1 / 2 = c x l- 1 ' 2 fj,K, (3) 

with cx given by c sync = 2.1, c/ ree = 13.7, and Cdust = 13.5 at 100 GHz. 
In addition, one has chi-u — 8-5 and chi-c — 20.6 for the Dust+Free-free 
components, uncorrelated and correlated (respectively) with HI. 

Of course, by their very construction, these normalisations are only appro- 
priate for intermediate and small scales, about i ^ 10. The corresponding 
fluctuation levels at 3.0 mm and 1.4 mm are shown in figure 7. By using these 
normalisations and our spectral model, we can compute the rms per 7 degree 
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beam at every frequency and check that the H I-correlated component indeed 
provides a good fit with the results of the |Kogut et al. 1 1996(1 analysis, as is 
shown in figure 3. 



3 Extragalactic foregrounds 



Many detailed studies have been devoted to effects that may blur the primor- 
dial signature of the CMB anisotropies. Gravitational lensing by mass con- 
centrations along the light ray paths may for instance alter the detailed map 
patterns and add a stochastic component. Or photons passing through a fast 
evolving potential well might be redshifted (Rees-Sciama effect). Secondary 
fluctuations might be generated during a reionisation phase of the Universe. 
For all this processes, the answer plees and Sciama, 1968| , |Ostriker and VisE~- 
aiac, 198"E| , |Seljak, 1996| , |Jafle and Kamionkowski, 1998] | is that their impact 
is quite small at scales corresponding to £ ~ 1/9 ^ 1000, and can easily be 
accounted for at the analysis stage of CMB maps. Further fluctuations may 
also be imprinted, in case of a strongly inhomogeneous reionisation; see e.g. 
Aghanim et al. |1996|| for further discussion. We shall ignore these effects in 



the following. 



3.1 Infrared and radio sources background 



Here we attempt to evaluate what might be the contribution from the Far In- 
frared Background produced by the line of sight accumulation of extragalactic 
sources, when seen at the much higher resolution foreseen for future CMB 
experiments. 

Given the steepness of rest-frame galaxy spectra longward of ~ 100 /im, pre- 
dictions in this wavelength range are rather sensitive to the assumed high- z 
history. Indeed, as can be seen from figure 4. a, a starburst galaxy at z = 5 
might be more luminous than its z = 0.5 counterpart because the redshifting 
of the spectrum can bring more power at a given observing frequency than 
the cosmological dimming, at least in the A ^ 800 /im range [e.g. Blain ancl 
Longair, 1993|1 (precisely the range explored by the HFI instrument aboard 
Planck). This "negative K-correction" means that i) predictions might be 
relatively sensitive to the redshift history of galaxy formation ii) variations 
of the observing frequency should imply a partial decorrelation of the galaxy 
contribution, i.e. one cannot stricto sensu describe this contribution by spa- 
tial properties {i.e. a spatial template) and a frequency spectrum (but see 
below). This is compounded by the fact that this part of the spectrum is not 
well known observationally, even at z ~ 0, and very little is known about the 
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redshift distribution of faint infrared sources. 




Fig. 4. Effect of redshifting a template starburst galaxy spectrum (left). Longward 
of ~ 300 /im, the flux of such galaxies are similar at all z ^ 0.5. 

Estimates of the contribution of radio-sources and infrared galaxies to the 
anisotropies of the microwave sky have till recently relied on extrapolations 
from redshift z = all the way to a (large) assumed galaxy formation onset 
redshift. Given the uncertainties described above, this makes the reliability of 
the predictions of this type of modelling hard to assess in the v > 100 GHz 
range. On the other hand, predictions of galaxy formation models in the UV 
and optical bands have received a lot of attention in the last few years; they 
rely on substantially more involved semi-analytical models of galaxy forma- 
tion which provide a physical basis to the redshift history of galaxies ||Whitej 
In short, one starts from a matter power spectrum (typi- 



and Frenk, 1991 



cally a standard, COBE-normalised CDM), and estimates the number of dark 
matter halos as a function of their mass at any redshift (e.g. using the Press- 
Schechter approach). Standard cooling rates are used to estimate the amount 
of baryonic material that forms stars. The stellar energy release is obtained 
from a library of stellar evolutionary tracks and spectra; computations may be 
then performed to compare in detail with observations. This approach, despite 
some difficulties, has been rather successful, and many new observations fit 
naturally in this framework. 

This type of physical modelling has recently been extended to the far-infrared 
range by puiderdoni et al. |1998|| . Given the stellar energy release, |Guiderdom 



et al. 1 1998|| estimated the fraction reradiated in the microwave range using a 
simple geometrical model, which yields an infrared luminosity function at all 
redshifts. This, together with an assignment of synthetic spectra to a given in- 
frared luminosity, allows predictions of the numbers and fluxes of faint galaxies 
at any frequency. To constrain the fraction of heavily reddened objects as a 
function of redshift, they used the z ~ IRAS data at 60 /zm [[Lonsdale et al 
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1990] to normalise their model, and they further assumed that the isotropic 
background discovered by Puget et al. 1 1996(1 is indeed the long-sought Cosmo- 
logical Infrared Background (CIBR) due to the accumulated light of infrared 



sources. 



Guiderdoniet al. [1997J selected a small set of possible redshift histories satisfy- 
ing this (integral) constraint, and predicted that ISO observations at 175 /im 
could help breaking the remaining degeneracy in the model. Such observa- 
tions were very recently done fKawara and al., 1997| , |Clements et al., 1997] , 
Dole et al., 1999| and the number of sources found is best described by one 



model of this family (their model E) corresponding to a fairly large fraction of 
very obscured objects. As shown by figure 5. a, the ISO Hubble deep field data 
at 15 /im ||01iver et al., 1997j1 also agree with the predictions of this model. 
Even the first SCUBA determination ||Smail et al., 1997] , plain et al., 1998] 
at 850 /im seems to be reasonably fitted by model E. In short, this model 
is successful in predicting the latest source counts over a broad range of fre- 
quencies. This gives us confidence that the contribution of infrared galaxies 
to CMB measurements performed in the same frequency range (that of the 
HFI, v > 100 GHz) can be assessed with decent accuracy. The situation is 
different for the redshift distribution of the sources which is more sensitive 
than the number counts to the details of the evolutionary scenarios. The cur- 
rent observational constraints are so far limited to the spectroscopic follow-up 
of IRAS sources at 60 /im which essentially highlights the IR properties of 
galaxies in the local universe. It is only when redshift follow-up of these new 
far-infrared catalogs will have been completed that new generations of models 
can be developed and refined assessment of the infrared source contamination 
can be made. 



In what follows, we shall assume that resolved sources have been removed, and 
we focus on the spatial and spectral properties of the remaining unresolved 
background. That can be done only by assuming a specific instrumental con- 
figuration giving the level of detector noise and a geometrical beam for each 
channel. Indeed the variance of the shot noise contributed by N sources of flux 
S randomly distributed on the sky is a 2 = NS and the corresponding power 
spectrum is C(£) = N/ (4tt)S 2 . For a population of sources of flux distribution 
per solid angle dN/dS/dQ the power spectrum is, after removing all sources 
brighter than S c 



C{t) = ^= T ^S 2 dS. (4) 

^beam •'O CLzLCLo 

1 /2 

If the cut S c is defined as S c = q (a 2 oni + aj ust + a 2 nstvum + a 2 CMB ...) with 
q fixed to 5, this formula allows to derive iteratively the confusion limit o" con f- 
This is the standard deviations of the fluctuations of the unresolved galaxy 
background, once all sources with S > 5a tot (i.e. including the other sources of 
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fluctuations, from cirrus, ddust, detector noise oWnim, the CMB, <Jq UB ) have 
been removed. 
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Theoretical estimates from model E of Guiderdoni et al. 1997| . (1) HFI wave- 



bands central frequency in GHz. (2) Beam full width half maximum in arcmin. (3) 
1 a instrumental noise for 14 month nominal mission. (4) 1 a fluctuations due to 
cirrus at Nhi = 1-3 10 20 cm~ 2 (level of the cleanest 10 % of the sky). The fluc- 
tuations have been estimated following Gautier et al. (1992) with P(k) oc fc~ 2 ' 9 
and Po.ioo^m = 1-4 10- 12 ^ 100Aim . (5) la CMB fluctuations for AT/T = 10" 5 . 
(6) la confusion limit due to FIR sources in beam f2 = 0p WHM , defined by 
G conf = (Jo hm S 2 (dN/dS)dS Q) 1 / 2 . The values a con f and Sn m = qa to t have been 
estimated iteratively with q = 5. (7) a to t = (c 2 ns + a conf + a cir + ^cmb) 1 ^ 2 - Here 
a c i r is for Nhi = 1-3 10 20 cm -2 . (8) Surface density of FIR sources for Sn m = 5atot- 



Table 1 summarises the E-model results for the HFI instrument Guiderdoni 



1998 



The main feature is the very high level of the confusion limit, which 
translates in modest number of detected sources. One should realize though 
that this estimate is very pessimistic since it assumes a very naive source 
removal scheme by a simple thresholding. In practice one would at least use 
a compensated filter (e.g. by doing aperture photometry) which removes the 
contributions to the variance of all long wavelength^ A much larger number 
of sources would then be removed since the counts are very steep (as a matter 
of example, |Guiderdoni et al. |1998|| predict with the same model more than 
120 000 sources at 350 /im (860GHz) at flux > 100 mJy, a detection threshold 
suggested by the recent work of |Hobson et al. |1998|| ). One should thus regard 
the confusion level derived above as an upper limit which very likely exceeds 
the real residual unresolved background from sources. 



3 A compensated filter (whose integral is naught) essentially nulls the contribution 
of fast decreasing power spectra contributors, like dust, whose contribution to the 
variance is concentrated at large scales. And it also decreases the amount of white 
noise from the detector and the unresolved sources themselves, since then sources 
have to stand out above a reduced threshold. 



16 



The spectral dependence of the confusion limit can be approximately modelled 
as that of a modified black body with an emissivity oc z/ ' 70 and a temperature 
of 13.8 K. A more precise fit (at the 20% level) yields at v > 100 GHz for the 
equivalent temperature fluctuation power spectrum (using relative bandwidths 
of 0.25) is given by 



£C(£) 



1/2 



7.1 10~ 9 

s ie/2.53 _ I 



0.16 N sinh x 



x 1 



x 



0.3 



[K], 



(5) 



with x = hu/2kT = z//(113.6GHz). This gives a value of C{£) 1 ' 2 = 0.005 /iK 
at 100GHz. At lower frequencies [y < 100 GHz), this model yields instead 



£C(£) 



1/2 



6.3 10~ 9 (0.8 - 2.5x + 3.38a; : 



sinh z 



x 



X 1 



m 



(6) 



Correction Footor (xy)/((x 2 }<y'»'« 
500x500 1.5'x1.5'E model, Poisson trunc 5.00 < 



1( - 









X (/im) 




is ; 




60 : 




x 175 : 










', IRAS, ISO & SCUBA 





-1 -2 

lo gio s [ J y] 




199^ 



Fig. 5. a) Number counts prediction from the E model of Guiderdoni et al. 
the available data (see text), b) Cross-correlation coefficients between maps based 
on the E model. 

Figure 4 shows that the relative fluxes in different observing bands can vary 
depending on the redshift of the object. Conversely, different bands do not 
weight equally different redshift intervals. The previous analysis thus does not 
tell us whether the fluctuation pattern at a given frequency is well correlated 
with the pattern at another frequency. Of course, the answer will depend on 
the resolution of the maps and the noise level. To answer this question for the 
[AAO 1 1998|| , maps of the galaxy contribution were generated at 30 different 



frequencies, as in |Guiderdoni et al. |1996|| and [Hivon et al. [1997]] , but for model 



E. As shown in figure 5.b, the cross-correlation coefficients of the maps (once 
the 5a t ot sources have been removed) is better than 0.95 in the 100-350 GHz 
range, better than 0.75 in the 350-850GHz, and it is still ^ 0.60 in the full 
100-1000 GHz range. This indicates that one can treat the IR background from 
unresolved sources reasonably well as just another template to be extracted 
from the data, with a fairly well defined spectral behaviour^, at least in the 



Note though that the model may not do full justice to the diversity of spectral 
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range probed by the HFI. 



While the model above might well be the best available for bounding the 
infrared sources properties, it does not take into account the contribution 
from low-frequency point sources like blazars, radio-sources, etc . . . which are 
important at the lower frequencies probed by MAP and the LFI . Using only 
that model would under-estimate the source contribution at low frequencies. 
On the other hand, |Toffolatti et al. |1998|| made detailed predictions in the 
LFI range which give 

tCltfl* 5.7sinh>/113.6) 

V '> ~ (y/l. 5) (4.75-0.185 log(v/1.5)) ^b I I 

This yields C(£) 1/2 =0.02 /iK at 100GHz. We shall also use that prediction 
when analysing the expected performances of MAP0. Here again, we assume 
that this unresolved background has well-defined spectral properties. 



In order to obtain a pessimistic estimate of the total contribution from sources, 
we shall consider in the following two uncorrelated backgrounds from sources, 
the one contributed by IR sources, as described by the equations (5) & (6), 
and the one contributed by radio- sources, as described by equation(7). The 
corresponding levels are compared with the other sources of fluctuations in 
figure 7 and they are quite low (at least at £ ^ 1000) as compared to the 
expected fluctuations from the CMB or from the Galaxy. 



3.2 Clusters of galaxies 



Hot ionised gas along the line of sight generate distortions of the CMB fluctu- 
ations when the incoming photons scatter off electrons and get a shift in fre- 
quency. This translates in an apparent temperature decrement in the low fre- 
quency (Rayleigh- Jeans) side of the spectrum, and a temperature excess in the 
Wien side. This is called the Sunyaev-Zeldovich (hereafter SZ) effect ||Zeldovichl 



and Sunyaev, 1969[| . The intensity variation of the CMB, AI U /I U = y x f(x) 



is controlled by the Compton parameter 

ko~T 



m P c 2 



T e (l)n e (l)dl, 



shapes of the contributing IR galaxies; if this turns out to be the case, we could for 
instance treat this contribution as an additional noise contribution to the detectors, 
as was conservatively assumed in the Planck Phase A study. 
5 Although in that case the remaining unresolved background should be somewhat 
higher due to the lower sensitivity and angular resolution of MAP. 
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where T e and n e stand for the electron temperature and density. The spectral 
form factor in the non relativistic limit depends only on the adimensionnal 
frequency x = hv /kT C MB according to 



/(*) 



xe 



x 



e x + 1 



4 



(9) 



During the Planck phase A study, |Aghanim et al. 1996|| devised a model 



for generating maps of the Sunyaev-Zeldovich effect (both thermal and ki- 
netic) and analysed the capabilities of Planck in detecting clusters of galax- 
ies [[Aghanim et al., 1997a| , see their Figure 1.6]. Since then, the model was 
improved by using better estimates of the counts (and of the peculiar veloci- 
ties), and generalised to encompass various cosmological models. The counts 
were derived from the Press-Schechter mass function 



Press and Schechter. 



1974 ], normalised using the X-ray temperature distribution function derived 



from |Henry and Arnaud 1 199 If data as in |Viana and Liddle |1996|| . The follow- 
ing power spectrum of the fluctuations due to the Sunyaev-Zeldovich thermal 
(hereafter SZ) effect from clusters of galaxies was evaluated by analysing maps 
of the Compton parameter y generated with these counts. 
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Fig. 6. Power spectra of the Sunyaev-Zeldovich thermal effect for different cosmo- 
logical models: standard CDM (solid line), open CDM (dashed line) and lambda 
CDM (dotted line) are shown. 

For the standard CDM model, it was found that the y fluctuations spectrum 
is well-fitted in the range 20 < £ < 4000 by 

£(£+l)C e = a ysz £[l + b ysz £}, (10) 

with a ysz = 4.3 x 10~ 15 and b ysz = 8.4 x 10~ 4 . Of course, the values of the 
fitting parameters in equation(lO) depend on the assumed cosmological sce- 
nario. For an open model (f2 = 0.3), one finds instead a ysz = 4.6 x 10 -15 and 
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b ysz = 2.2 x 10~ 3 . These are quite small differences, as illustrated by figure 6. 
This is natural since most of the contributions come from relatively low red- 
shifts while the counts have precisely been normalised via z ~ observations. 
Equation(lO) for standard CDM yields 



'C 1 / 2 = 0.27 



'1 + 8.4 x 10" 



V 2 , X 

fiK (11) 



for the temperature fluctuation spectrum at 100 GHz. 

On small angular scales (large I), the power spectrum of the SZ thermal effect 
exhibits the characteristic dependence of the white noise. This arises because 
at these scales the dominant signal comes from the point-like unresolved clus- 
ters. On large scales (small /), the contribution to the power comes from the 
superposition of a background of unresolved structures and extended struc- 
tures. The transition between the two regimes occurs for I ~ l/b sz that is when 
the angular scale is close to the pixel size of the simulated maps. One should 
note though that this approach neglects the effect of the incoherent superposi- 
tion of lower density ionised regions like filaments. While the contributions in 
the linear regime are easy to compute ||Persi et al., 1995| , the non-linear con- 



tributions must be evaluated using numerical techniques IPersi et al, 1995} 



Fogosyan et al., 1996||. While the modelisation above is quite sufficient for 



the purposes of this paper, detailed simulations of the component separation 
should benefit from simulated maps of the SZ effect with more realistic low 
contrast patterns. 

The kinetic SZ effect due to the Doppler effect from clusters in motion along 
the line of sight (hereafter l.o.s) cannot, of course, be spectrally distinguished 
from the primary fluctuations. This additional contribution is about an order 



of magnitude smaller than the thermal SZ effect ||Aghanim et al., 1997a| ,p|. 
We shall thus neglect this contribution to the temperature fluctuations in the 
following analysis. 



4 Comparing contributions to the microwave sky 



4-1 Description of the detectors noise 



In the course of the measuring process, detector noise is added to the total 
signal after the sky fluctuations have been observed, which can be described 
by a convolution with the beam pattern. In order to directly compare the as- 
trophysical fluctuations with those coming from the detectors, it is convenient 
to derive a fictitious noise field "on the sky" which, once convolved with the 
beam pattern and pixelised, will be equivalent to the real one fKnox, 1995 
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Modelling the angular response of channel i, Wi(9), as a Gaussian of FWHM 
9i, and assuming a white noise detector noise power spectrum, the sky or "un- 
smoothed" noise spectrum is then simply the ratio of the constant white noise 
level by the square of the spherical harmonic transform of the beam profile 



cm 



exp 



(l+l) 2 ' 



02Q2 



; exp( 



2^2 In 2 



:i2i 



with 



2 sin 



9; 



2 v / 8hi2 / 



and c noise 



affti = af x 2vr[l - cos(^/2)], if a { stands for the 1-a AT/T 
sensitivity per field of view. The values of c no i se for different experiments of 
interest can be found in Table 2. 



While convenient, the noise description above is quite oversimplified. Indeed 
detector noise is neither white, nor isotropic. The level of the noise depends 
in particular on the total integration time per sky pixel, which is unlikely to 
be evenly distributed for realistic observing strategies. In addition, below a 
(technology-dependent) "knee" frequency, the white noise part of the detector 
spectrum becomes dominated by a steeper component (typically in 1/f, if 
the frequency / is the Fourier conjugate of time). This extra power at long 
wavelength will typically translate in "stripping" of the noise maps, a common 
disease. Of course, redundancies can be used to lessen this problem flWright 
or even reduce it to negligible levels if the knee frequency is small 



1996 



enough ||Janssen and al., 1996|, pelabrouille, 1998 
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Table 2 

Summary of experimental characteristics used for comparing experiments. Central 
band frequencies, u, are in Gigahertz, the FWHM angular sizes, Ofwhm, are in 
arcminute, and AT sensitivities are in fiK per Ofwhm x Qfwhm square pixels; 
the implied noise spectrum normalisation c no i se = AT^fwhm) 1 ^ 2 , is expressed in 
/xK.deg. 



21 



4-2 Angular Scale Dependence of the Fluctuations 

Figure 7 compares at 30, 100, 217 and 857 GHz the power spectrum of the 
expected primary anisotropies (in a standard CDM model) with the power 
contributed by the galactic emission (eq. [3]), the unresolved background of 
radio and infrared sources (eq. [5, 6, 7]), the S-Z contribution from clusters 
(eq. [10]), and on-sky noise levels (eq. [12]) corresponding to the PLANCK 
mission. It is interesting to note that even at 100 GHz, the dust contribution 
might be stronger than the one coming from the synchrotron emission, at 
least for levels typical of the best half of the sky. Since point processes have 
flat ("white noise") spectra, their logarithmic contribution to the variance, 
£{£ + 1)C{£) oc £ 2 , increases and becomes dominant at very small scales. 

4-3 Frequency Scale Dependence of the Fluctuations 

Given the power and frequency spectra of the microwave sky model, one can 
then compare the rms contributions per beam at any frequency for any 
experiment by 

a(v l ) = J2W+l)C e (v l )w l (£) 2 , (13) 

l 

where W{ stands for the transform of the beam profile at the frequency V{. 
Figure 8 gives an example for the Planck case. 

4-4 Angular- Frequency Dependence of the Fluctuations 

An extensive series of angular power spectra at many frequencies (such as 
those in figure 7), or of amplitudes of fluctuations as a function of frequency 
for a series of angular scales (such as figure 2), provides a detailed picture of 
the behaviour of the various components. Still, it is also illuminating to retain 
the full dependence of the fluctuation amplitudes on frequency and angular 
scale to build a more synthetic view of the fluctuations "landscape" . 

Figure 9 shows the contours in the v — £ plane where the fluctuations, as 
estimated from £(£+1)C(£)/2tt, reach one tenth (left) or one hundredth (right) 
of the large scale COBE level [£{£ + l)C(£)/2n = (30/iK) 2 /10 ~ (10/iK) 2 ]. 
These contours map the three dimensional topography of the fluctuations of 
individual components in the v — £ plane. The synchrotron component (when 
expressed in equivalent temperature fluctuations) defines a valley which opens 
towards large £, since £ 2 C(£) decreases with £ as l/£. The H I-uncorrelated 
component defines a shallower and gentler valley (see also fig.2.b), while the 
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Fig. 7. Contributions to the fluctuations of the various components, as a function 
of angular scale, at different frequencies. The thick solid black line corresponds to 
the CMB fluctuations of a COBE-normalised CDM model. The dots, dashes and 
dot-dashes (in red, blue, and green) refer respectively to the HI correlated, HI uncor- 
related, and synchrotron emissions of the galaxy. The light blue triple dots-dashes 
displays the contributions from unresolved radio sources, while the purple line cor- 
responds to the unresolved contribution from infrared sources. Long dashes show 
"on-sky" noise level (see text, eq. (12)). On the 30 GHz plot, the highest noise level 
corresponds to MAP, and the lower one to the LFI. On the 100 GHz plot, the 
highest noise level corresponds to the LFI, and the lower one to the HFI. In the 
two remaining plots, the noise level corresponds to the HFI case. 
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Expected contributions in the model 




10 100 1000 

v [GHz] 

Fig. 8. Standard deviation per beam versus frequency for all the relevant model 
components in the Planck case. 

H I-correlated emission creates a high frequency cliff. The large I end of the 
valley is barred by the point sources "dam". In any case, even in the case of 
Planck, the noise leveQ is higher than the contribution from the unresolved 
sources background. 

The heavy black line shows the path followed by a stream lying at the bottom 
of the valley, i.e. it traces the lowest level of total fluctuations (but noise). 
Its location confirms that v ~ 100 GHz is the best frequency for low-£ mea- 
surements. At £ 200, the optimal frequency moves to higher values and is 
determined by the minimum of the fluctuations from unresolved sources and 
clusters; indeed, the dotted line shows the minimum location if clusters' fluctu- 
ations are ignored. Taking them into account shows that the optimal frequency 

6 The noise levels plotted were computed using a linear interpolation between the 
columns of the summary table 2. 
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Contours levels of l.C(l) 1/s = lOfsK Contours levels of l.C(l) 1/2 = 'suK 
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Fig. 9. Contour levels of the different components of our model in the angular 
scale-frequency plane. The levels indicate when the power spectra of the various 
components reach (10/iK) 2 and (3/iK) 2 - respectively one tenth and one hundredth 
of the level of the COBE detection; the filled side of the lines shows where the 
corresponding process is higher than the level line (with the same colour coding 
than in previous plots, but for the HI correlated and uncorrelated components which 
have been replaced here by the total dust and free- free emissions). Thus the central 
area is the region where a CMB signal at 1/10 or 1/100 of the COBE level will 
dominate over all other components. The central dashed line materialises the locus 
of weakest overall contamination by foregrounds. As before the galactic components 
corresponds to the median sky level. The noise level shown is for Planck. 

is £ 200 GHz for high resolution work around i ~ 1000, rather close to the 
zero of the Sunyaev-Zeldovich effect at 217 GHz. Of course the exact opti- 
mal value depends on the relative weight of the unresolved backgrounds from 
sources and that from clusters, and cannot be very precise at this stage. Still, 
the optimal range is ~ 150 — 200 GHz, which is within the range accessible 
to the most sensitive detectors available today (typically bolometers can span 
the 100 to 1000 GHz range). In the case of Planck, one can expect that 
the most stringent constraints on the high-£ part of the CMB spectrum will 
come from the 143 & 217 GHz channel, not the lower frequency ones of lower 
angular resolution. 

In order to better identify the dominant contributions in a given experiment, 
Figure 10 shows the regions of the v — i plane where the contribution from 
a given process tC v {€) accounts for half of the total including the noise, 
I J2pC p (£) (i.e. of the contribution per logarithmic interval centred around 
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Figures available with full version at 

ftp '■ 1 1 ftp.iap.fr /pub/ ' fromjusers /bouchet/ 'wienerl '.ps.qz 

Fig. 10. Each hatched areas delimits the regions where a given component accounts 
for 50% of the total contribution to the variance including noise, in a particular 
experiment. Each panel has been restricted to the frequency range that the experi- 
ment span. Inner contours show when a component accounts for 90% of the signal, 
and the central contour (horizontal hatches) delimits the region where the CMB 
anisotropies should contribute 99% of the signal, a) on the left is shown the MAP 
case, b) the right panel corresponds to Planck. 

that scale I to the variance of the measurement at frequency v). The inner 
contours show when that process accounts for 90% of the total (and 99% for 
the CMB). This figure thus shows the areas in the v-l plane where a given 
component will be dominating the measurements. 

Two cases were considered for illustration. On the left, the noise level and 
the frequency range correspond to the MAP experiment which uses pas- 
sively cooled HEMT detectors, while the Planck case which combines ac- 
tively cooled HEMTs and bolometers is shown on the right. In the latter case, 
reaching 10% accuracy on the CMB anisotropies over all modes between a few 
and a thousand, in the cleanest 50% of the sky, should be easily accessible. A 
1% accuracy can only be obtained though in a substantially smaller angular 
range, if no foreground subtraction is performed. We now turn to that issue. 



5 Separation of components 



Our goal in this section is to introduce a simple tool, which we denote as 
the "quality factor" of an experiment, that will allow in §6 a comparison of 
experimental set-ups, once the effect of foregrounds is included. We start by 
establishing notations for describing our model of the sky (§5.1) and derive 
in §5.2 a formulation of the separation problem. General properties of linear 
inversion methods are described in §5.3. The Wiener filtering method is de- 
scribed in §5.4, and §5.5 gives the properties of this method, while §5.6 shows 
that a simple indicator allows capturing most aspects of the inversion. Simple 
examples and more technical points are discussed in the final sub-section. 
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5.1 Physical model 



We note T{y, e) the flux at a frequency v in the direction on the sphere 
referenced by the unit vector e. We make the hypothesis that T{y, e) is a 
linear superposition of the contributions from N c components, each of which 
can be factorised as a spatial template T p {e) at a reference frequency (e.g. at 
100 GHz), times a spectral coefficient g p (v) 

F(v,e) = y Eg p {v)T p (e). (14) 
P =i 

This factorisation is not stricto sensu applicable to the background due to 
unresolved point sources, since at different frequencies different redshift ranges 
dominate, and thus different sources. But even at the Planck-HFI resolution, 
this is, as we showed earlier, a weak effect which can be safely ignored for our 
present purposes. 

Let us also note that the factorisation assumption above does not restrict the 
analysis to components with a spatially constant spectral behaviour. Indeed, 
these variations are expected to be small, and can thus be linearised too. For 
instance, in the case of a varying spectral index whose contribution can be 
modelled as u a ^ e '7 p (e), we would decompose it as u a [l + (a(e) — a) \nv]T p (e). 
We would thus have two spatial templates to recover, T p (e) and (ot(e) — a)T p (e) 
with different spectral behaviours, oc v a and oc v a \nv respectively. But given 
the low expected level of the high latitude synchrotron emission, this is unlikely 
to be necessary. This simple trick may still be of some use though to describe 
complex dust properties in regions with larger optical depth than assumed 
here (see also the alternative way proposed by |iegmark~ 1998|| ). 



5.2 The separation problem 



During an experiment, the microwave sky is scanned by arrays of detectors 
sensitive to various frequencies and with various optical responses and noise 
properties. Their response is transformed in many ways by the ensuing elec- 
tronics chains, by interactions between components of the experiments, the 
transmission to the ground, etc. The resulting sets of time-ordered data (TOD) 
needs then to be heavily massaged to produce well-calibrated pixelised maps 
at a number of frequencies, each with an effective resolution, and well specified 
noise properties. We assumed here that all this work has already been done 
(for further details, see |Janssen and Gulkis 1 1993^1 , |Wright |1996|| , [legmark 
Delabrouille 



1997 



1998H, Delabrouille et al. 119981 



We thus focus on the next step, i.e. the joint analysis of noisy pixelised sky 
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maps, or signal maps Sj, at N v different frequencies in order to extract infor- 
mation on the different underlying physical components. It is convenient to 
represent a pixelised map as a vector containing the N p map pixels values, the 
j-th pixel corresponding to the direction on the sky e^. Let us denote by t p the 
unknown pixelised maps of the physical templates T p . We thus have to find 
the best estimates of the templates t p (ej), which are related to the observed 
signal by the N v x N p equations 



e i) = H / dv Vi(v - Vi) g p (v) J2 Wi(ej - e k ) t p (e k ) + ni(ej) (15) 
v =\ Ju=0 k=i 

= J2i v i' k 9p(v)] [wi*t p (e)) +ni(ej), (16) 
P =i 



where * stands for the convolution operation. This system is to be solved 
assuming we know for each signal map its effective spectral and optical trans- 
mission, Vi[y — vi) and Wi(e) and that the covariance matrix Nj = (ninj^ of 
the (unknown) noise component rij has been properly estimated during the 
first round of analysis. 



Since the searched for templates t p are convolved with the beam responses, 
Wi, it is in fact more convenient to formulate the problem in terms of spheri- 
cal harmonics transforms (denoted by an over-brace), in which case the con- 
volutions reduce to products. For each mode £ = {£, m}, we can arrange 
the data concerning the N u channels as a complex vector of observations, 
y(£) = { si (£), s 2 (£), • • • , Sn„(£)}, and we define similarly a complex noise 
vector, h(i) = {ni (£) , r rv^{l) , . . . ,n/v^(£)}, while the corresponding template 
data will be arranged in a complex vector x(£) of length iV c . We thus have to 
solve for x, for each mode t independently, the matrix equation 



y(£)=A(£)x(£)+b(£), (17) 

with A ip {£) = [vi * g p {v)\ Wi(£). Note though the independence of the £ modes 
would be lost in case one uses other functions that spherical harmonics (e.g. 
to insure orthogonality of the basis functions when all the sky is not available 
for the analysis fGorski, 1994j| ). Then instead of solving N p systems of N v 
equations, one will have to tract a much larger system of N p x N v equations 
(possibly in real space rather than in the transformed space). For simplicity 
we treat each mode separately in the following. 
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5.3 Linear Inversions 



Any linear recovery procedure may be written as an N c x N v matrix W(£) 
which applied to the y(£) vector yields an estimate x.(£) of the x(£) vector, 
i.e. 

x(£) = W(£)y(£). (18) 

In the following, we show that the matrix R = WA plays in important role in 
determining the properties of the estimated templates by the selected method, 
whatever it might be. Indeed, since x = W [Ax + b], we have 

(x) = R (x) and (xx'\ = R (xx f ) R f + WBW f , (19) 



where x^ is the Hermitian conjugate of x (i.e. x^ is the complex conjugate of 
the transpose, x T , of the vector x), and B = (bb^ stands for the (suppos- 
edly known) noise covariance matrix. We denote by (Y) the average over an 
ensemble of statistical realisations of the quantity Y. The covariance matrix 
of the reconstruction errors, E, is given by 



E = ((x - x)(x - x) 1 ") = C - RC - CR 1 " + (xx 1 ") , (20) 
where C = (xxM stands for the covariance matrix of the underlying templates. 



An estimate of this covariance matrix C should of course be estimated from 
the data itself. Since the elements of the covariance matrix of the recovered 
templates (assuming they have zero mean) may be written as 

./■„.,■: ) =R PP (x p x* q ) iC + ££ R pr (x r x*) R\ q + [WBW+l , (21) 



V q j 

rjtp s^q 



pq 



we define the intermediate variable x p = x p /R pp in order to write 

(x p S%) = (x p x*) + B, (22) 

where B is then an additive bias. This shows that an unbiased estimate of the 
covariance matrix of the unknown templates, C, may be obtained by 

m =+e. 

C = (xxt) = — — y xx 1 " - B. (23) 

m=—l 

This expression is only formal though, since B contains some of the unknowns. 
Still, this form allows to generalise the calculation of [[Knox, 1995|| of the 



expected error (covariance) of the power spectra estimates derived in such a 
way from the maps. Indeed, 
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Cov 



a 



pi 



a 



pq 



a 



pq 



2£ + 



pq 



-a 



pq 



(24) 
(25) 



Once this expression is developed, there are many terms with fourth moments 
of the form (x r (£, m) x s (£, m)* x t {£, ml) x u {£, m')). 

If we assume that all the x p are Gaussian variables, then we can express these 
fourth moments as products of second moments, and recontract the sums to 
finally obtain 



Cov 



a 



pq 



— {C pq + B pq ) 2 + 2£_\_l ( PI + Bpi) 

— 2Bpq((C pq + Bpq) + B 2 q — Cp q 

2 ,„ „ ,o 2 



2£+ 1 



{Cpq + Bpg) 2 



Xp Xq 



2£+l \RppRq-q, 



(26) 
(27) 

(28) 



This expression shows the extra contribution of noise and foregrounds to the 
cosmic variance. It can be re-expressed and simplified once a specific inversion 
method has been specified. 



5.4 Wiener filtering 



One particular choice of W may be obtained by minimising the figure of merit 



X 



: (€) = (y-Ax)t IT 1 (y-Ax), 



(29) 



for each mode i. The advantage of this simple inversion by x 2 minimisation is 
that it requires no prior information on the unknown x, only some knowledge 
on the statistical properties of the noise maps through their covariance ma- 
trix B. An assessment of this simple method will be given in a sequel paper 
(hereafter paper II) concerned with actual numerical simulations of the sky 
and its analysis. 

Another possibility for linearly separating the components is to apply "Wiener" 
filtering techniques since [Bouchet et al. |1996b|| and [legmark and Efstathiou 



1996|| showed how to generalise this well-known procedure (see e.g. Punn et al 
Bond 



1995|| for an application to COBE-DMR data, and Zaroubi et al 



1995|| for galaxy surveys data) to the case of multi-frequency, multi-resolution, 



data. In the following, we rederive in a self-contained way this Wiener filter- 
ing approach, with a number of extensions which will prove useful for our 
purposes. 
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The inversion by x 2 minimisation is designed to best reproduce the observed 
data set by minimising the residuals between the observations and the model 
(eq. [29]). Another possibility is to require a statistically minimal error in 
the recovery, when the method is applied to an ensemble of data sets with 
identical statistical properties. One then obtains the "Wiener" matrix W, by 
demanding that the vector of differences between the underlying and recovered 
templates, e = x — x, be on average of minimal norm 

e 2 (£) = (|x - x| 2 ) = ((x - x) t (x - x)) (30) 

In practice, we require that the derivatives of 



5 2 (£) = (|x- W(Ax + b)| 2 ) (31) 
= (W pv A vp , - VX^VA/p" - VO + W pv W p»' ( b "K>) ( 32 ) 

versus all W pv (and their complex conjugate) be zero for each £ (with repeating 
indices implying summation in(32) and (33)), which yields 



W, 



pv 



A up , (xptx*,,) A v , v n + (b u bl,) = A vlp , (x p 'X p ) . 



(33) 



In matrix form, the wiener matrix is then given by 



W = CA T 



ACA T + B 



(34) 



if we denote by C the covariance matrix of the templates, C, 



vv 



j-.pj-.pi 



Alternatively, W f = ACA T + B AC (since A is real, A f = A T ). 

Since the Wiener matrix depends on C, i.e. on the covariance matrix of the 
templates, it is clear that this makes use of prior knowledge to construct opti- 
mal weights. One can thus think of Wiener filtering as a "polishing" stage of a 
first inversion (e.g. as in § 5.3) yielding an estimate C of the covariance matrix 
needed to have optimal weights. A practical implementation is described in 
paper II. 



5. 5 Characteristics of the Wiener Separation 



The definition (33) of Wiener filtering implies that 



xxM = (xxM = (xx^ 



(35) 
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As a result, the expression for the covariance matrix of the reconstruction 
errors E = (eeM reduces to the simple form 



E = [I - R] C, 



(36) 



where I stands for the identity matrix, and R = WA as before. Alternatively, 



E = C 



I-Rt 



The general expression of the covariance matrix of the templates given by 
eq. (19) may be simplified by using the specific Wiener matrix properties of 
eqs. (35) and (36), which yields 



xx 1 



t\ = R/xxt\. 



(37) 



Both equations(36) and(37) clearly show that the closer R is to the identity, 
the smaller are the reconstruction errors. 

The diagonal elements of E are of particular interest since they give the resid- 
ual errors on each recovered template, 



£p (1 Rpp)C^ 



pp- 



(38) 



To break down this overall error into its various contributions, we start from 
eq. (32) 

Ep — (1 — Rpp) 2 Cpp + Rpp'Rpp"Cplp" + WpyWpylByyl . (39) 

p'+p+v" 

By using eq. (38), we get 

1 



Bp — — — ^ Rpp'Rpp"Cp' P " + W pi ,Wpv'B uv i, 
ttpp p'jtp^p" 



(40) 



which describes the power leakage from the other processes and the co-added 
noises. 

In addition, the error on the deduced power spectra given by eq. (28) takes 
for Wiener filtering the simple form 



Cov 



C 



AC 



Rpr C r q 

~R R~^ 

1 ^PP 1 l qq 



2i+l 

which will prove particularly convenient to compare experiments. 



(41) 
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5.6 Case of uncorrelated templates 



In this theoretical analysis, we can always assume for simplicity that we have 
decomposed the sky flux into a superposition of emissions from uncorrelated 
templates, so that 

(xpx$) = 5 pp , C pp (£) = 5 ppl C p (£). (42) 

In our galactic model, it means in particular that we consider the templates for 
the H I-correlated and H I-uncorrelated components rather than those for the 
dust and free-free emissions. In actual implementations though it might prove 
more convenient to relax this assumption and consider correlated templates 
but with simpler spectral signatures (see Paper II). In the uncorrelated case, 
the previous expressions (37, 38 & 40, 41) may be written as 



x p \ 2 ) = Q p {\x p \ 2 ) = Q P C P (£) (43) 

e\ = (1 - Q P )C P = 1/Q P £ R\ v ,C pl + W W W VV ,B VV , (44) 

v'+v 



where Q p = K pp stands for the trace elements of WA. Thus Q p tells us 

(1) how the typical amplitude of the Wiener-estimated modes x are damped 
as compared to the real ones (eq. [43]). 

(2) the spectrum of the residual reconstruction error in every map (eq. [44]); 
note that this error may be further broken down in residuals from each 
component in every map by using the full R matrix. 

(3) the uncertainty added by the noise and the foreground removal (oc 1/Q P — 
1, eq. [45]) to cosmic (or sampling) variance which is given by C pp /V2lTT 
(this result only holds though under the simplifying assumption of Gaus- 
sianity of all the sky components). 



Given these interesting properties, we propose to use Q p as a "quality factor" 
to assess the ability of experimental set-ups to recover a given process p in 
the presence of other components; it assesses in particular how well the CMB 
itself can be disentangled from the foregrounds. 

Section 6 below is devoted to "practical" uses of the quality factor, in con- 
junction with our assumed sky model. 



5. 7 Discussion 
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5. 7. 1 Simple examples 



This "quality" indicator generalises the real space "Foreground Degradation 
Factor" introduced by Podelson and Stebbins |1994j| . It may be viewed as 



an extension of the usual window functions used to describe an experimental 
setup. This can be seen most easily by considering a noiseless thought ex- 
periment mapping directly the CMB anisotropics with a symmetrical beam 
profile w(8). Then the power spectrum of the map will be the real power 
spectrum times the square of the spherical harmonic transform of the beam, 
< x p (£) >= w(£) 2 < x p (£) >. The spherical transform of Qp^ 2 (£) is then the 
beam profile of a thought experiment directly measuring CMB anisotropies. 
The shape of Q p (£) thus gives us a direct insight on the real angular resolution 
of an experiment when foregrounds are taken into account. 

In addition, one often considers (e.g. for theoretical studies of the accuracy 
of parameter estimation from power spectra) a somewhat less idealised ex- 
periment which still maps directly the CMB, with a beam profile w(6), but 
including also detector noise, characterised by its power spectrum Cn- Let us 
suppose this is analysed by Wiener filtering. We have to solve y = Ax + n, 
with 

A T = {w,w,w,w, ...}, (xx^ = C x , B = C N x I. (46) 
The previous formulae then leacQ to 

Qx {i) = C x (£)+w-*C N (£)/N c (47) 
A4 = = ^7 ( C S) + ^Y 2 C N (£)/N C ) , (48) 

where N c stands as before for the number of channels (measurements). Equa- 
tion(47) tells us that the Wiener-estimated modes (as seen by the ratio of the 
estimated to real spectra) are damped by the ratio of the expected signal to 
signal + noise (the noise power spectrum being "on the sky" , cf. §4.1, with the 
noise power spectrum being the noise power spectrum per channel divided by 
the number of channels). In brief, what Wiener filtering does when the noise 
strength is getting larger than the signal is to progressively set the estimate 
of the signal to zero in an attempt to return only true features. And of course 
we recover the traditional error on the estimated power spectrum added by 
noise found by |Knox 1 19951 . Thus determining the quality factor allows to 



find the equivalent ideal experiment often considered by theorists; it provides 



7 By using [I + fil}' 1 = I - (3/(1 + (3N C ) 1, if 1 stands for the N c x N c matrix 
with all elements equal to 1, and J for the identity matrix. 
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a direct estimate of the final errors on the power spectrum determination from 
a foreground model and the summary table of performance of an experiment. 



5.7.2 Reduced variables 

Note that the expressions found previously could be further simplified by 
considering reduced templates, x' which are uncorrelated and of unit variance. 
Indeed, let us perform a Cholesky decomposition of C, i.e. C = LL T where L 
is a lower triangular (real) matrix. The reduced templates are then defined by 
x' = L _1 x = L T C x x (which implies, as desired that the reduced templates 
are "orthonormalised" , /x'x't\ = I/C" 1 <W\ C~ 1T L = I). 



Then the problem to solve may be written as 

y = Sx' + b, (49) 



where we have defined S = AL. Since (yy / = SS T + B, SS T is the useful 
signal covariance matrix (and S may be seen as it's Cholesky decomposition). 
In that case, 



W' = S f [SS f + Bl 1 (50) 



E' = ((x'-x')(x'-x') t ) = I-R' (51) 
R' = L- 1 RL = S f [SSt + BpS. (52) 

This form shows most clearly that the Wiener filter derived is indeed a gener- 
alisation of the usual signal to noise weighting (eq. (52)). On the other hand, 
this convenient form for theoretical analysis might not be so practical since C 
might be rather ill-determined from a first step of analysis. Thus paper II will 
rather use the formulae obtained in the previous sections. 



6 Comparing experimental set-ups 



6.1 Examples of Wiener matrixes 



Once the instrument, through the A matrix, and the covariance matrix C 
of the templates are known, the Wiener filter is entirely determined through 
equation (34). Figure 11 offers a graphical presentation of the resulting values 
of the Wiener matrix coefficients of the CMB component when we use our sky 
model (assuming negligible errors in designing the filter). It shows how the 
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different frequency channels are weighted at different angular scales and thus 
how the v — £ information gathered by the experiment is used. 

In the MAP case, most weight is given to the 90 GHz channel. It is the 
only one to gather CMB information at i £ 600, and it's own weight becomes 
negligible at i ^ 1000. In the LFI case, the main CMB channel is the 100 GHz 
one, which gathers information till £ ~ 1500, while the 70 GHz channel is of 
some help till £ ~ 1100. The two other channels only contribute at £ ^ 500. 
The situation is different in the HFI case, since no single channel dominates 
at all scales. The 143 GHz channel is dominant till £ ~ 1400 and useful till 
£ ~ 2000, while the 217 GHz channel becomes dominant at i ^ 1400 and 
gather CMB information till £ ~ 2300. The 100 GHz channel contribution to 
the CMB determination is only modest, and peaks at £ ~ 800. 

Once the information from the LFI and the HFI are considered jointly, the 
previous HFI weights are barely affected, and the LFI weights appear to 
be sizeable only in the 100 GHz channel, albeit at a fairly reduced level as 
compared to the 100 GHz channel of the HFI (which is itself only a modest 
contributor to the CMB determination). As a result, this already suggests that 
the impact of the LFI on the high £ HFI measurements (£ > 200) will be in 
controlling systematics and possibly in better determining the foregrounds. Of 
course, one should note that this analysis assumes the foregrounds are known 
well enough to design the optimal filter. Thus the LFI impact might be greater 
in a more realistic analysis when the filter is designed with the only help of 
the measurements themselves. 



6.2 Effective windows and beams 

Figure 12 gives the effective £-space window of the experiment for each com- 
ponent, Q p {£), to be compared with the individual windows of each channel, 
Wi(£). For the CMB, it shows the gain obtained by combining channels through 
Wiener filtering. For MAP, the effective beam is nearly equal to it's 90 GHz 
beam, while the improved noise level of the LFI allows it to do much better 
than it's "best" beam would suggest (in particular at £ ~ 1000). None of these 
two experiments can separate the SZ or the infrared sources background con- 
tribution. As expected the Galactic foregrounds are poorly recovered except 
on the largest scales (small £) where their signal is strongest (we assumed 
spectra oc £~ 3 ). The only exception is that of the H I-correlated component 
which is of course very well traced by the HFI. When the LFI and the HFI 
are combined, the main improvement appears to be on the SZ contribution 
thanks to the greater spectral level arm offered to distinguish it from other 
foregrounds. 
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Fig. 11. CMB Wiener matrix elements at t > 200 for MAP, the LFI, the HFI and 
the full Planck mission. Each bin in £ is the average over one hundred contiguous 
values. The frequency bins are centred at the average frequency of each channel, 
but their width in not representative of their spectral width. 

Figure 13 shows the inverse spherical harmonics transforms of the Q p {£), i.e. 
the effective angular beams. They show the effective point spread function 
by which the underlying maps of the component emissions would have been 
convolved once the analysis is complete. One can see in particular that the 
FWHM of the CMB beam is ~ 16.5' for MAP and ~ 7.5' for Planck (and 
even better for the dust components). 

6.3 Maps reconstruction errors 

We now estimate the spectrum of the reconstruction errors in the map. Fig- 
ure 14 compares these residual errors per individual I mode with the typical 
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MAP-X experiment (1998) LFI - PLANCK 




lll'l - PLA.MCa PLA.\CK cx :.:M ir---:nt 




Fig. 12. Square of the £-space effective windows (i.e. Q p ) for the LFI, HFI, and 
PlanckAs usual, black is for the CMB, red, blue, and green are for the Galactic 
components, yellow is for the SZ contribution. The transform of the channel optical 
beams are shown for comparison as dotted orange lines. 

amplitude of the true signal. The first thing to note is that the error spectrum 
cannot be accurately modelled as a simple scale-independent white noise (as 
expected from the v — I shape of the Wiener matrix for the CMB) . We also see 
that the LFI improves on MAP mostly by decreasing the noise contribution, 
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Fig. 13. Corresponding real space effective beams (negative parts are denoted by 
dashes). As usual, black is for the CMB, red, blue, and green are for the Galactic 
components, yellow is for the SZ contribution. 

but at the largest scales where the galactic contributions are also lowered. The 
LFI taken jointly with the HFI helps in reducing the \cw-i galactic emissions, 
and the power leakage from the radio-sources background. The contributions 
per logarithmic bin of i to the variance are compared in figure 15. It is in- 
teresting to note that in the Planck case the largest contribution but noise 
comes from SZ clusters in the range 50 £ £ £ 500, to be then superseded by 
the leakage from the radio-sources background. 

6.4 Foreground and noise induced errors on the CMB power spectrum 

Finally, we can estimate the uncertainty added by the noise and foreground 
removal to the CMB power spectrum. Figure 16 shows the envelope of the l-a 
error expected from the the current design of MAP (green), the LFI (blue), 
and the HFI or the full Planck (red), as well as a "boloBall" experiment 
meant to show what bolometers on balloon might achieve soon; the experi- 
mental characteristics used in this comparison may be found in table 2. Note 
that there has been no "band averaging" in this plotQ, which means that 
there is still cosmological signal to be extracted from the HFI at i ~ 2500. 
The message from this plot is thus excellent news since it tells us that even 
accounting for foregrounds, Planck will be able to probe the very weak tail 
of the power spectrum and allow breaking the near degeneracy between some 
of the cosmological parameters. 

8 The band averaging would reduce the error bars on the smoothed C(£) approxi- 
mately by the square root of the number of multipoles in each band, if the modes 
are indeed independent. 
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Fig. 14. CMB reconstruction error contributed by each component (with the same 
line coding than above) and their total in black. The integral of l{l + l)e 2 /2tt 
would give the reconstruction error of the map. The dashed-dotted line shows for 
comparison C(£) 1 / 2 for the CMB. For Planck, the total of the reconstruction error 
is nearly constant at all I and corresponds to e(l) ^ 0.025 /xK, a factor greater than 
six below the MAP case. 
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Fig. 15. CMB reconstruction errors compared with the input spectrum in the MAP 
(left) and Planck case (right). 

The second panel of the figure shows the variations induced by changing the 
target CMB theory. In this flat universe with a (normalised) cosmolog- 

ical constant of 0.7 was assumed. What MAP can say on the third peak would 
then be much reduced, while on the contrary the HFI would probe accurately 
one more peak. This illustrates the fact that the errors on the signal will of 
course depend on the amplitude of the signal; it is a reminder of the fact that 
the numbers given so far are meant to be illustrative. As already stressed, 
the merit of this approach is to offer a quantitative tool for comparisons. The 
impact of changing the foreground model will be addressed in section §6.7. 



6.5 Wiener versus "naive" error estimates: 



In the previous section we estimated the error bars reachable on the CMB 
power spectrum assuming a sky model and a separation method. In effect we 
found an equivalent ideal experiment described by an "on-sky" noise power 
spectrum [w~ 2 C^} e ff = CqmbWIQcmb — !]• Numerous papers have been de- 
voted to estimating the possible accuracy of determinations of the cosmological 
parameters from CMB experiments. To forecast this accuracy, various recipes 
were used to derive an estimate of this effective on-sky noise spectrum from 
tables of performance of the experiment such as table 2. One option has been 
to assume that all frequency channels but one are used to remove the fore- 
grounds contribution from the "best" CMB channel, i.e. the one closest to 100 
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GHz, or the most sensitive one. A less radical proposition by Pond et al. [1997 



has been to assume that all channels but Nk eep are used to clean the Nk eep 
closest to 100 GHz, which are then assumed to contain only the CMB signal 
and noise, with an effective "on-sky" power spectrum of this noise [w~ 2 C^] e ff 
given by 



[w 2 C N } eff fr[ Wi C N>i 

This is a generalisation of the formula (47) to unequal noise power spectra, 
Cj\r,i in different channels. This approach does not tell how many and which 
iVfc ee p channels can be considered clean once the iV c — Nk eep others have been 
somehow used for foreground removal. 

Figure 17 compares the fractional error ACcmb/Ccmb expected from Wiener 
filtering and different prescriptions, namely using the best channel or all chan- 
nels in some frequency range combined according to equation (53). As expected 
we see that the error in the "naive" estimates of the uncertainty of the CMB 
power spectrum goes to zero at small £, when cosmic variance dominates over 
the noise, and at large I when detector noise dominates all foreground contri- 
butions. If only the best channel (dashes) of the experiment is retained (i.e. 
v = 90,143 GHz respectively for MAP & Planck), the error is of course 
always over estimated, i.e. optimal use of all channels by Wiener filtering does 
much better than this naive estimator would suggest. 

The error in estimating the noise level is much decreased if the channels be- 
tween 55 and 300 GHz (i.e. the 60 & 90 GHz channels of MAP and the 70, 
100, 143 and 217 GHZ channels of Planck) are retained in the sum (53). 
A better estimate for MAP obtains if one also includes the 40 GHz channel, 
since one now underestimates the error by only the few per cent over the full i 
range. For Planck on the contrary retaining all channels between 55 and 300 
GHz already underestimates the errors by 10 % at i ~ 1500. Assuming that 
all Planck channels between 35 and 400 GHz can be cleaned with the help 
of the remaining ones further increases the error on the CMB uncertainty to 
more than 15 % at i ~ 1800. 

More importantly, this error is scale-dependent, i.e. it shows that no "naive" 
prescription gives a fair estimate of the CMB power spectrum uncertainties 
at all scales. This is not surprising, given the shapes of the Wiener matrixes 
we saw in figure 11. These shapes indeed illustrate that an optimal estimate 
of the CMB component obtains by weighting differently as a function of scale 
the informations coming from different channels, according to their angular 
resolution, and to the levels of their detector noise and of their foregrounds 
contaminations. 

The results above show that the relative error between Wiener based and other 
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error estimates of the CMB spectrum are relatively small, but they are scale 
dependent, and the same naive prescription might either underestimate or 
overestimate the errors depending on the experiment considered. This invites 
exercising some restraint in using results based on these simple prescriptions, 
like the usual comparison of accuracy of cosmological parameters determina- 
tions by different experiments. Of course, our Wiener based estimates cure 
only some of the limitations of standard analyses since many simplifying (and 
unrealistic) assumptions remain. The main limitations of course arises from the 
idealised description of the measurement process (eq. (17)) and the assumed 
perfect knowledge of the response matrix, A, and the correlation matrix, C. 

MAP-X experiment. (1998) PLANCK experiment 




500 1000 1500 2000 1000 2000 3000 4000 5000 

1 1 

Fig. 17. Ratio versus scale of naive error estimates of the CMB power spectrum to 
the one derived from the quality factor. This ratio is given both for MAP (left) and 
the full Planck (right). 

6.6 Experiment optimisation and robustness 

One may wonder whether the CMB residuals will vary linearly with global 
variations of the noise figures of an experiment, as could be produced for 
instance by variations of the operating temperature, or length of operation of 
the experiment (allowing different number of sky coverages), or simply overly 
optimistic or pessimistic estimates of the detectors performances. Figure 18 
shows the variation of the rms of the total residual contribution from noise 
and unsubtracted foregrounds to a CMB map as a function of the noise level 
compared to the nominal one. The square of this rms value is obtained by 
integrating the residual spectrum of figure 14 till a maximum £ max . Varying 
d-max gives a concise description of the contribution of various scales to the 
map error. 

In the nominal case, a/a = 1, one sees that for MAP and Planck respec- 
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MAP-X experiment (1998) PLANCK experiment 

Fig. 18. Performance variations for global variations across channels of the detector 
noise level (e.g. with the number of sky coverages done). The yellow, green, and blue 
lines correspond to restricting the integration to a maximum £ ma x of respectively 
500, 1000, 1500. The red line corresponds to no restriction at all. 

tively the residuals are 21 and 2.5 /xK for "low" -resolution maps retaining only 
the modes till £ max = 500. For medium resolution maps with £ < 1000, these 
errors increase to 41 and 5.5 fiK, and high-resolution maps with £ < 1500 
have residuals of 50 and 9 fiK respectively. If all £ are retained, the rms of the 
residual are then 52 and 13 fiK. The difference between MAP and Planck 
in the increase of the residual with the number of modes retained reflects 
the difference in angular resolution of the maps produced by Wiener filtering 
(figures 12 and 13). 

For "low" and medium resolution maps, the residuals for Planck vary nearly 
linearly with the noise level. This is the high-£ part which is most affected by 
variations of the noise level. In the MAP case, the stronger curvature of the 
lines indicate that it would be more affected by sub-nominal performance of 
the detectors, and conversely that it would benefit more of an increase of the 
mission duration; a decrease of the noise level by a factor of two (e.g. 4 years 
instead of one year of operation) would for instance reduce the residuals for 
£ < 1000 from 41 to ~ 30 yK. 

The quality factor approach may also be used to test whether the failure of 
a single channel of an experiment would be critical for the outcome of the 
experiment, i.e. whether it is robust. Here we estimate the variation of the 
CMB residuals in different "crippled" configurations for Planck. Figure 19 
shows that i) the residuals are dominated by noise ii) there is nearly no im- 
pact at all if a single Planck channel is lost except for the 143 & 217 GHz 
channels hi) loosing the HFI would more than double the residuals, but the 
LFI completed by the HFI 217 or 353 GHz channel would already do rather 
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Fig. 19. a) CMB residuals (in a map with all modes till £ ma x = 1500) retained) in 
various configurations for Planck. The black curve displays the total, and the other 
lines detail the various contributions from noise (orange curve) and foregrounds 
with the same colour coding than figures 7 or 8 (i.e. red, blue, and green curves 
correspond to the Galactic residuals, yellow is for the thermal SZ, and the light blue 
and purple curves are for the unresolved sources backgrounds), b) Total residual in 
different configurations, for different £ m ax (as in figure 18, the yellow, green, and blue 
lines correspond to restricting the integration to a maximum £ max of respectively 
500, 1000, 1500, while the red line corresponds to no restriction at all), c) Table of 
configuration indexes in the adjacent figures. 
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6. 7 Uncertainties of the sky model 



Figure 16 illustrated the fact that the deduced 1-cr errors on the CMB spec- 
trum depend on the amplitude of that signal by comparing the results for 
a standard and a "Lambda-CDM" model. Given the observed variations one 
may worry that all the results obtained so far might critically depend on the 
details of the sky model assumed. This is potentially frightening since clearly 
the situation regarding foregrounds emission might be quite more complicated 
than assumed here (e.g. some free-free may locally be emitted at higher tem- 
perature, the backgrounds from unresolved sources might be higher, etc). 

We have thus considered variations on the basic sky model used so far. For 
the three galactic components, we considered cases when each power spectrum 
might level off and become oc £~ 2 at i > 200, resulting in an increase of power 
by a factor of 10 at £ = 2000 as compared to our basic assumption of a constant 
slope of -3 at all scales. We also considered that the level of the free-free like 
emission uncorrelated with the Hl-emission might be four times lower. We 
also increased independently the normalisation of each unresolved foreground 
component by a factor of four. We allowed each of these 6 variation to occur 
in combination with all the others, thereby creating 64 different foreground 
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models. 



Figure 20 shows the resulting 100 GHz power spectra in 5 of these 64 config- 
urations, as well as the corresponding quality factors for MAP and Planck. 
As anticipated, the effective transmission of a foreground component is im- 
proved when and where the amplitude of that component is increased. But 
the transmission of the other components does not appear affected. This is 
not really surprising either, since Wiener filtering is a signal to (signal+noise) 
weighting, where (signal+ noise) obtains for the contribution from detectors 
noise and of of all the astrophysical components. The lack of sensitivity of 
the transmission of all components but the one varied simply reflects the fact 
that the dominant terms are the CMB and detector noise, which are precisely 
those which were kept untouched here. 



MAP-X experiment (1998) 



PLANCK experiment 




1000 1500 2000 2500 

1 




Fig. 21. Effect of variations in the sky model (see text) on the deduced 1 — a error 
area for the CMB power spectrum. Both in the MAP (left) and Planck (right) 
case, the 64 variations investigated barely affect the result because the dominant 
contributions to the signal in the channels of most weight for the CMB determination 
are the CMB itself and the detector noise. 

Figure 21 displays the extent of 1 — a error for the CMB power spectrum 
in the case of MAP and Planck for the 64 sky models explored. Nearly 
all curves are superimposed, which demonstrates that the deduced errors are 
rather robust to variations in the details of the sky model. This is simply 
because the dominant contributions to the signal in the channel of most weight 
for the CMB determination remained the CMB itself and the detector noise. 
Thus if the main features of the sky model proposed here are correct, it is 
unlikely that most results obtained in this paper will be much modified. 
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7 Conclusions 



In order to assess the expected level of foreground contamination to be ex- 
pected in CMB experiments we have developed a simple sky model based on 
current knowledge. In this model, the microwave emission of the sky is the 
superposition of the CMB, three galactic emission mechanism, in addition to 
two unresolved backgrounds from sources, one associated with radio-sources 
and the uncorrelated other one arising from infrared sources. Further fluctu- 
ations are imprinted by the Sunyaev-Zeldovich effect of clusters. While not 
known precisely, we believe that the non-CMB components are reasonably 
well constrained by current data supplemented with modest theoretical ex- 
trapolations. 

We then analysed aspects of the data analysis of multi-frequency multi-resolution 
observations. Wiener filtering is shown to allow building optimal linear esti- 
mates of the spatial templates of the physical processes, as defined by their 
spectral behaviour. We computed specific examples of Wiener filters for our 
sky model and several experimental set-up. The shape of the CMB Wiener 
filter in the v — i space allows to build intuition on how the information from 
different frequency channels is combined at different angular scales to min- 
imise the residuals in a CMB map. We further showed as a function of scale 
the fraction of residuals which arise from detector noise and from partly un- 
subtracted foregrounds. Integrating these residuals show for instance that the 
rms reconstruction error from modes I < 1000 should be ~ 40 fiK for MAP 
and smaller than 6 fiK for PLANCK. 

We also introduced a simple measure of the effective resolution of multi- 
frequency multi-resolution experiments in the presence of foregrounds. It's 
harmonic transform is the effective point spread function by which the un- 
derlying maps of the component emissions would have been convolved once 
the analysis is complete. This effective transmission or window function is 
directly related to the expected uncertainty on the power spectra determina- 
tions once foregrounds are accounted for (under the assumption of Gaussian 
foregrounds though). Typically, MAP will probe well the CMB power spec- 
trum till i ~ 900 (1000) and Planck till I ~ 2300 (2500) for a standard 
(or a Lambda) CDM model. We find that simpler estimates make a scale- 
dependent error in estimating the spectrum uncertainty and that depending 
on the experiment they might yield an over- or an under-estimates. 

As further examples of the usefulness of this simple analytical framework, we 
computed how the rms reconstruction error is affected by global noise level 
variations in an experiment, as induced for instance by a variation in the 
mapping duration of the experiment. MAP is relatively sensitive to such vari- 
ations, while the variation is closer to linear for Planck. We also analysed the 
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robustness of Planck to various failures leading to the loss of some channels. 
The most important channels are the 143 and 217 GHz, but even their total 
loss would result in less than a 10 % increase in the residual rms. Loosing all 
HFI channels but the 217 GHz one would only yield a 20 % increase, but it's 
loss would then nearly double the residual rms! Clearly, optimisation trade-off 
should not result in decreasing the sensitivity of the 217 GHz and neighbouring 
channels. 

Of course, assessing the effect of foregrounds relies to some extent on the 
realism of the foregrounds model. We considered 64 sky models differing by 
what appears to be the least well established features. We found that the 
quality of the CMB recovery should be basically unaffected, mainly because 
the dominant high-£ contributions in the channels of most weight for the CMB 
determination are the CMB itself and detector noise. This excellent news 
should nevertheless be taken with a grain of salt. Indeed, it was always assumed 
a nearly optimal filter could be built, which entails some knowledge of the 
foregrounds (i.e. their very existence, their power spectra, and their frequency 
dependence). Paper II will address some of these points by means of direct 
numerical simulations. We can already note though that a strong foreground 
at least in some channels should be no problem since the needed information 
will be easily obtained from the data itself. Forgetting a very weak foreground 
in all channels should not lead to great errors. The difficult case will of course 
be to make sure that no unsuspected foreground is sizeable mostly in the 
primary CMB channels. 

But much new data is becoming available to improve on our current knowledge 
of the microwave sky, either from the analysis of data from the ISO satellite, 
of ongoing H a and radio surveys, and of many ground and balloon experi- 
ment. And of course the future generations of CMB satellite themselves were 
designed to much increase this knowledge by mapping most of the sky in a 
large frequency range and with much improved noise levels as compared to 
COBE. The analysis above thus makes us rather confident that foregrounds 
should not prevent the future CMB satellites from finally unveiling the CMB 
power spectrum with high accuracy. 
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